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Introduction 


A method for calculating three dimensional separated flow around 
bodies similar to ground vehicles was investigated in collaboration 
with the Office National d' Etudes et de Recherches A6rospatiales in 
Chatillon. This investigation was a continuation of research started 
in 1975 in which we have analyzed many results obtained in wind and 
water tunnels and proposed [15] a two dimensional approach based on a 
wake-vortex representation. 

In the present article the investigation has been extended to 
three dimensions, and a nonstationary, inviscid scheme based on the 
method developed by Rehbach [1] is proposed. 

The obstacle is represented by singularities and the wake is 
modeled \’/ith the aid of vortex particles emitted from the separation 
line, which is assumed to be known in this case. 

The calculations carried out for axisymmetric bodies show the 
evolution of the vortex sheet's structure as a function of the angle 
of the flat base area. Wind tunnel visualizations confirmed the 
results. The sheet, -which is of a separated two dimensional nature 
for an upright base area, takes on a split sheet structure for 
greatly inclined base areas [8] , with lateral rollup present. 
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The numerical results presented constitute a first approximation 
of the problem. They especially shov7 the necessity of an in-depth 
study of three dimensional separation. 

Chapter 1 — The Physical Problem /l. 

1.1 Introduction 


A physical analysis of the base area flow past simple bodies 
related to ground vehicles was conducted in the Lelarge-Renault wind 
tunnel at the Institut A6rotechnique at Saint-Cyr-1 'Ecole . The 
experiments' main goal was to complete the recent information 
obtained by other authors. (Morel, Maull [2] in 1976) 

The body investigated (figure 1.1) was composed of a common 
forward part made up of a Rankine half -body [3] with a diameter 
0 = 160 mm and length L = 700 mm. Fitted to it are flat base areas 
with angles a = 0® , 30®, and 60®. The assembly is mounted at zero 
incidence and sideslip with aid of streamlined supports. Height H 
can take on two values; H/J2f = 2 and E/0 = 0.65. 
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The analysis was based on; 

— streamline visualizations using paste 
— wake visualizations with bubbles (EMIBUL technique) 

— tracings of boundary layer profiles at the base area edge 

[ 6 ]. 


The results presented here are those for H/0 = 2. 

The boundary layer readings make possible an approximate 
calculation of the intensity of the wake vortex emitted at the 
separation line 2)’- 



Fiaure 1 ,2 
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The circulation 6r on contour C, which is defined in figure 
1.2 is, during the period of time fit: 

-V • ndo (1*1) 

If one supposes that the flow is locally two dimensional, 
one obtains, with: 



and 


with 


curlV.n 



d/ => udt 



/ 


1 dr r\j 
V2dC“ 




( 1 . 2 ) 


where Ug is the velocity outside the boundary layer. It has been 
assumed here that the velocity u^ at the base area is negligible. 

1.2 Experimental Results 

The observations and photographs reveal the following elements: 
(Rl ^ 1.2 X 10^) . 

1 . 2 . 1 a = 0 


The vorticity discharge emitted by the boundary layer is 
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constant along the separation line *5^ (figure 1.3 and table 1.1, in 
which the values of Ug/Voo, 6, and K given by equation (1.3 are noted). 


Table 1 . 1 


9 

Ue/Voo 

^(mm) 

K 

0° 

1.03 

10 

0.53 

30° 

1.02 

10.5 

0.52 



Figure 1.3 

The vorticity vector u associated with each vortex ring is 
parallel to the yz plane, and the vortex sheet I is symmetric about 
the x-axis. Finally, the streamlines ij) remain parallel to the x-axis 
and are perpendicular to line S. The lift F^ of such a body when 
it is located far from the ground is zero and its drag is situated at 
point a of the curve Cx = f (a) (traced in figure 1.8, which is 
excerpted from [2]). 

1.2.2 a = 30° 


The streamline visualizations show that the streamlines ^ 
are perpendicular to the line ^ . The slight base area inclination 
causes dissymmetry in the wake, but its two dimmensional structure 
is conserved. Representative drag point b in figure 1.8 confirms 
this. 

1.2.3 g = 60° 

i) The structure of the vortex sheet is greatly modified even 
though, as shown in photo 1, the streamlines remain perpendicular to 
S. (See figure 1.4.) 
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on secondaire 


e 


c 


Key: a) dividing point 

b) main vortex 

c) secondary vortex 


d) region swept by main vortex 

e) accumulation line 
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streamline visualizations of the base area reveal secondary 
rollups similar to those observed in swept wings (Werle [7]). These 
rollups are illustrated in figure 1.5, which is based on photographs. 


/ 1. 4 



Figure 1.5 


Key; a) wake rollup 

These observations were recently confirmed by H. Werle through 
visualizations made in a water tunnel on a body similar to the one 
described in figure 1.1. 

In addition, three”dimensional effects are preponderant. The 
vortex sheet splits at A in figure 1.5, and a rolled up wake is 
emitted all along contour . 

In order to satisfy this scheme, it is necessary that the 
vorticity vector u emitted be tangent at B to line §| and evolve 
continually from B towards A in such a way that it is practically 
parallel to at A. 

ii) The photographic data obtained with the bubble technique 
(see photos 2 and 3) clearly show the formation of the two lateral 
vortices. The flow remains separated in the median plane of the 
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base area, probably with the formation of a bulb at A, as the down- 
stream point B behaves like the following edge of a wing. 



PHOTO 2 photo 3 


iii) This last remark is confirmed by the readings of static 
pressure on the wall at the median line (figure 1.6). The high pres- 
sure in the downstream region shows that the extrados streamline 
is strongly deflected in the direction of V by the intrados flow. 
(See also figure 1.9, which was taken from [2].) 



figure 1.6 














(■Cpb) 



figure 1.9 


excerpted from [2}, T. MOREL 


Chapter 2; Theoretical Basis 
2.1 The Problem 


It is proposed to solve the following problem, representing the 
flow of an incompressible, inviscid fluid, defined by the equations 
[4] : 


||-Vxv^V = -grad(£+iV^) 

(2.1) 

O 

II 

> 

(2.2) 

and the boundary conditions: 


V.?|3=o 

(2.3) 

V— Voo , -* oo 

(2.4) 

^ in the equations designates the velocity vector, 

p is the 


pressure, and ^ the mass density. The slip condition (2.3) is 
applied to solid obstacles inside an unbounded domain with uniform 


/ 2.1 


10 




velocity at infinity. 


By taking the curl of (2.1) so as to intorduce the vector 
w = curl the preceding system can be written in the purely 
kinematic form: 


where 


designates the particle derivative in Helmholtz’s equation (2.5). 

2.2 Calculation of the Velocity Field /2 . 2 

In the problem (2.5), (2.6) with boundary conditions (2.3) and 
(2.4), the velocity field V can be represented by Green's function 
{9, 1]: 


(uI,gracl)V 

dl/_* 

db'vV. 0 


^=1^ f(y.graoj^ 


(2.5) 


( 2 . 6 ) 


uj^X3raol[l)dVq + |^divVej.gracl(l)dVG^ 


IT 


(2.7) 


This function gives an expression for V at any point P of a 
volume V inside a closed surface S whose normal n points inward. 


When point P belongs to volvime v, the second member of (2.7) 
has the value 4ir^p. 

It is equal to zero when P is outside volume v. 

Let there be two concentric surfaces Sq and Si (figure 2.1). 
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Let Vi be the volxome between Sq and and V 2 be the volume inside 
Si. 



2.2.1 /2.3 

Equation (2.7) written for a point P inside volume v^^ is; 

- |ffxv)*^(i)dire| - 

- |(nxV)x^(-p)d(T^ - |p.Y).Cjr^|ijd(rq = 4-TxVp (2.8) 


2 . 2.2 


since the surface Sj^q is that of a sphere whose radius R -»■ », 
equation (2.7) together with condition (2.4) yields for a point within 

SlO- 

-|(nxV)xp(i)clira, - 1 (n.v) 9raS(i)d(ro=+n\i, 

Sio 510 
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and relation (2.8) becomes; 


2.2.3 



(2.9) 


It is interesting to find any possible discontinuities in velocity 
across S^ 2 • ^o do this, (2.7) is written for P outside volume V 2 : 


1 4ivV2.3radJi|d->iGi 
- l"xV>grad[i]40ffl-|n.V>3rddW^^ 

^ 4 ?. ■' 'n ' ( 2 - 1 °) 


S21 

If it is assumed that, in volume V 2 ; 

div^= 0 (2.11) 

0 

Adding (2.9) to (2.10) v;hile including the boundary condition 
(2.2) results in; 


/2.4 


'hrVp= Woo + 



(a) 

(b) (2.12). 

(c) 


where (+) and (-) have replaced subscripts (12) and (21) . 


The integral equation (2.12) can be considered as the basis of 
the method of singular points. The problem posed by (2.5) and (2.6) 
and the boundary conditions (2.3) and (2.4) are replaced by (2.5), 
(2.12), and condition (2.3) on S. 
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2.2.4 Method of Singular Points 


i) Integral (a) of equation (2.12) expresses the contribution 

at P of all vorticity vectors u contained in voliime It should 

be noted that the volume integral only covers the regions of v^ 
where 5 does not equal zero. This makes its calculation easy when 
the domain v is unlimited. 

In addition [1] , as the fluid is as sinned to be incompressible, 
the calculations are carried out for the vector; 

( 2 . 13 ) 

where w represents the average vortex intensity associated with the 
particle of constant volume STf . 

ii) Integral (b) can be written, with = -n” = n; 



The local scalar quantity q = ft* - ’^”) is equal to the increase 
in the normal component of velocity across the surface element dag 
and therefore has the dimensions of a rate per unit area. 

The contribution of velocity corresponding to integral (b) is then 

'WVp.| = I - • 3rad[ij (245) 

When point P approaches point Q (figure 2.2), the contribution 
of the integral at the singularity is 2irqpft''' on the n"*" side and 
2 ivqpft“ = -27rqpft’^ on the n~ side, with the increase ft*(^^- ^“) equal 
to the local intensity q, 

iii) Integral (c) of (2.12) can also be written; 
and the local vectoral quantity y = n x (^ - ^“) is equal 


(2.16) 
to the 


/2.5 
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surface vortex density associated with the discontinuity in tangential 
velocity across surface S. 


— + 



The contribution of velocity corresponding to integral (c) is: 

iwV^=, ( 2 . 17 ) 

This is the classic Biot-Savart equation [10]. The integral ' s 
contribution at the singularity is 2Tr^xn'*’ on the n"*" side and 2n^ xtH” - 
-2Trfx^I+ on the n~ side, with the increase in the tangential component 
of velocity equal to t^xti. 


2.2.5 Application to Physical Problems 

The integral equation (2.12) is written: 



i) The hypothesis of continuity of tangential velocity on S, or /2 . C 
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leads to the "source" solution. It is known that this solution is 
particularly well adapted to external zero-lift Neumann problems. In 
contrast, the distribution of sources giving an imposed distribution 
of normal velocity on an internal surface S (internal Neumann 
problem), is generally not unique. (Hunt [5]). 

ii) The hypothesis of continuity for normal velocity on S, 

mm mm 

n ‘V = n 'V , leads to the "vortex" solution. 

This solution is suitable for problems with positive lift and 
makes it possible to take into account vortex sheets for incompressible 
flows of transient viscosity. The vortex intensities in the sheet are 
then directly derived from the intensities y on the surface of the 
body (Rehbach [1]). 

Let a volume element [11] 

JVq 3 = 

and w be the vorticity vector at a point Q or a surface S (figure 
2.3). The quantity e can represent the local thickness of the 
boundary layer and 5 the transported viscosity vector. 



The application of Guass's theorem to inf inetesiraal cylindrical 
volume element dv yields the following approximate result: 



(2.19) 


CU 5cl(r =s^n1xV^-f- n"xV^clo+ ^lateral surface 

U)£d(T« lateral surface 

Note that ^ =■ 

By making e^O (Reynolds ->■ ”) in such a way that the product 
ue remains constant, equation (2.19) yields: 

7 == ( 2 . 20 ) 

which is the vectorial quantity already obtained in paragraph 2.23 
with Green's function. 

Finally, (figure 2.4) if surface S is a closed surface on which 
the restriction 6j2f/5n = 0 (^'^•n^ = 0) everywhere is imposed^ and if 
the problem is resolved with the vortex solution, equation (2.19) 
shows that Sj?'*'/6n = 0. It can be deduced from this that the potential 
within the interior volume is constant and that ^“ = 0 everywhere [5] . 



The local intensity y = therefore is the image of the dis- 

continuity of tangential velocity on surface S when the Reynolds 
number approaches infinity. 


In these conditions the intensities 0 of the vortex particles 
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which we used to model the wake Z '{figure 2.4) and are followed with 
the aid of Helmholtz's Equation (2.5) can be determined directly with 
the aid of the vortex intensities ^ on surface S along the separation 
line . 

iii) In contrast, it is not necessary to include the idea of 
an infinite Reynolds number for the source solution. 

The vorticity discharge into the wake can in principle be cal- 
culated on the basis of physical data valid for a fixed Reynolds 
number. 


Three-dimensional emission models do not seem to exist at 
present. Approximate schsmes based on local two-dimensional approaches 
are used, and the vorticity discharge has the form (1.2) ; 



( 2 . 21 ) 


In equation (2.21) Ug and u^ represent the velocities on either 
side of 2) (figure 2.4); k is uniform at a given distance, see [12] 
and [1] . The main difficulty in the use of this model is in the 
choice of the unit vector e, whose direction is closely linked to the 
extent of the flow's three-dimensional character at the base area. 


Thus, for an upright base area, is everywhere tangent to the 
separation line, whereas for a highly inclined base area, the vector 
^ varies continuously along 2). • 


This point is numerically analyzed in the following chapter. 
2.3 Calculation of Pressure 


The pressure field is related to the velocity field by equation 
(2.1) , which, in the regions where curl ^ = 0, i.e. outside the 
vortex sheets, is written: 

C(t) (2.22) 


/ 2.8 
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where = grad 0 and C(t) is a constant function of time t. 

It will be assumed in the following discussion that a steady 
state on surface S has been achieved and that the potential 0 tends 
toward a finite limit, which is physically likely. 

2.3.1 Upright Base Area 

Let there be a cylindrical body whose axis is directed along the 
vector and which possesses an upright base area. 

Physical observations (Werle [7], Siriex [16]) show that down- 
stream from the base area there is a vortical recirculation region 
R2 (figure 2.5a) separated from region Rl by streamline C. 



figure 2.5 


Point R is a free reattachment point. 

It is known in addition that phy&ically and for Reynolds numbers 
on the order of 10^, the velocity at the base area is small. The 
pressure coefficients, reduced to p.oo and [Vco] , are on the order of 
-0.2 at the axis. This value remains practically constant over the 
base area, but depends on the configuration of the edge. 


/2.9 
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The constants and C 2 in regions Rl and R2 are therefore 
different/ and for points 1 and 2 located in the two regions: 


“ F~+ = fii 

2 

P 2 + ■= pi2 


If one examines a point 2 on the base area where V 2 = 0 and a /2 . 10 
point 1 in Rl where Vi = 0 also, one will therefore have: 

fLlEl =r 

k'- - id- 

or AKp =s • 

zC^ 

The quantity AKp is the irreversible energy loss due to viscous 
effects. This loss is not taken into account in a perfect fluid model 
leading to the result in figure 2.5b. It therefore seems difficult to 
obtain realistic base area pressure coefficients in this case. (See 
also [13] / p. 14. ) 

2.3.2 Inclined Base Area (a = 60”) 


In cases in which the sheet is split, the problem is different. 
For points on the base area which are not located near the separation 
line f there is a steady-state streamline in which (figure 2.6) : 


The pressure coefficients can then be calculated by: 

Kp= ■'-[ 



(2.23) 


(2.24) 


20 



/ 



Key: a) vortex wake rollup 

2.4 Calculation of Stress /2 . 11 


2.4.1 


When a velocity-pressure relation can be defined for all points 
in S, the aerodynamic force can of course by calculated by: 


F = I 


(2.25) 


2.4.2 


Referring to equation (2.23), equation 2.24 becomes: 



(2.26) 


In 

185-86) 


the vortex solution equation 
becomes: 


(2.26) (Milne-Thomson [11], 


PP 



(2.27) 


The velocity vector can be decomposed into: 


21 





=* Vco + + V/ojq 

where ^u>Q ^respectively represent the velocities induced in 

Q by the vortex layer ^ distributed on S and by the point singularities 
contained in the sheet 2 . 


The aerodynamic force I* is then: 

F, t Fi + Ft 

with 





(2.28) 


In ?Q and the conventional forms of induced lift and drag can 
be recognized. 


Chapter 3 ; Nvimerical Results 
3 . L Introduction 

The different calculation programs developed on the basis of the 
preceding theoretical elements have made it possible to numerically 
study three-dimensional separated flow around cylinders with flat 
base areas inclined at an angle a (figure 3.1) . 



Key: a) panel 


/3.1 
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The bodies are broken up into discrete quadrilateral panels. 
The control point for each of the panels is its center of gravity. 


The calculations presented here were carried out for bodies at 
zero incidence and sideslip located at a great distance from the 
ground, and with L/R = 3. 

The evolution of the vortex particles emitted at the separation 
contour D is calculated on the basis of Helmholtz's equation (2.5), 
which for nondimensional variables is written: 



where S = 1/V At is the Strouhal number associated with 1 and and 
the vortex release frequency n = 1/At. 

The choice of step At in the discrete equation associated with 
(3.1) determines the Strouhal nvimber of the numerical model. If the 
length 1 designates the average characteristic dimension of the panels 
located at [ sic ] , the value of S corresponding to a satisfactory 
evolution of the sheet is on the order of one. 

3.2 Numerical Results: Vortex Solution 


For this type of numerical model and for the results presented, 
the emission is made from a flap (figure 3.1a) located directly down- 
stream from Sb • 

The method for creating the vorticity vector tl is similar to 
the one proposed by Rehbach [1] for the study of flat wings of zero 
thickness. The adaptation to three dimensions is sketched out in 
figures 3.2a, b, and c. 


/3.2 
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figure 3.: 


For a = 0° and 30®, the sheet has a two-dimensional separated 
structure, and the vectors ^ associated with the particles emitted at 
the flap's following edge remain parallel to ^ (figure 3.4b). 

Finally, particularly for a = 0®, the velocities on the base 
area are very small, which is physically correct. 

For a = 60®, the emission lines traced on figure 3.5a show 
lateral rollups beginning to form. The flow is separated on the lower 
part of the base area and reattached at its upper part, in the median 
plane. 

It should be pointed out that with this model, the emitted 
vectors ^ remain practically parallel to the following edge of the 
flaps since the sheet is not split. The direct consequence of this 
is that the velocities obtained at the base area are insufficient, and 
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the pressure coefficients remain positive there. 


The improvements made recently by Billet [14] in the emission /3 
model (elimination of the flap, participation of the panels neigh- 
boring ) should make it possible to obtain good results, while 
at the same time probably causing the sheet to tear. 

3.3 Numerical Results: Source Solution 


For this type of numerical model, the emission is defined by 
equation (2.21), with the quantity u^ not being taken into account 
in the present shceme. The flap associated with the scheme described 
in paragraph 3.2 is suppressed, and the direction of the vector ^ is 


For a = 0® (figure 3.6) , the 
vector is parallel to [ sic ] . The 
result shows a zone of recirculation 
beginning to form downstream from 
the base area. The velocities on 
the base area, reduced to V», are on 
the order of 0.1, which is physically 
correct. 


arbitrary. 

3.3.1 
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3.3.2 


For a = 60®, several numerical tests are presented. Their goal 
is to analyze the role of the vectors' direction at emission, whose 
modulus is defined by equation (2.21): 

i) The vector ^ remains parallel to : The result in figure 

3.7, obtained after an 18-step calculation, confirms the one in 
figure 3.5. The sheet is not split; the velocities at the base area 
remain insufficient. 





Vecto: 

along 

18 -st( 

culat 





ficure 3.9 


ii) The vector ti varies continuously along contour ^ : The flow 

is two-dimensional in the downstream part of the base area (figure 3.8). 

The numerical results obtained are interesting, the sheet is split, and 
the lateral rollups are satisfactory. The velocities at the base area 
are greater, but insufficient at ^ . This is probably due to the fact 
that the split in the sheet is incomplete, with the emission line 
(marked with the symbol on figure 3.8) remaining in the median 
part of the base area. 

iii) Since the force exerted on the vector particle ^ has the 

form "P 'V ^ xV, the sheet is in equilibrium when ^ || For the prob- 

lem that interests us, the steady state was the only one desired, and 
it seemed useful to impose these conditions on the emission from the 
beginning of the calculation. 

The results obtained were similar to those analyzed in paragraph /3 . 8 
ii (figure 3.9) as the velocities at the base area were insufficient. 

It also appeared that the emission line (marked with the symbol 
"«"), which was located downstream from the base area, participated in 
the lateral rollup. This seems to be in contradiction with the exper- 
imental measurements, which show that this part of the base area behaves 
like the following edge of a wing (Chapter 2) . 

3.4 Conclusion 


The body of results presented here shows the possibility of 
calculating separated flow past thick bodies through the use of a 
theoretical model based on a discrete wake-vortex representation. 

The calculations carried out specify the profound changes in the 
wake when the base area angle is large. 

Our results constitute an as yet incomplete first approximation of 
the pressure on the base area and show that an in-depth analysis of the 
mechanism of three-dimensional separation is necessary. This work is 
in progress. 


30 



REFERENCES 


[1] Rehbach, C., "Ecoulements avec nappes tourbillonnaires" [Flow 
with Vortex Sheets] , La Recherche Aerospatiale , 5 (1977) . 

[2] Morel, T., D.J. Maull, and W.H. Hucho, Aerodynamic Drag Meehan - 
isms of Bluff Bodies and Road Vehicles , Plenum Press, 1978. 

[3] Milne-Thomson, L.M. Theoretical Hydrodynamics , Macmillan, 1938 
and 1950. 

[4] Lamb, Sir H. , Hydrodynamics , Dover Publications, 1945. 

[5] Hunt, B., "The Panel Method for Subsonic Aerodynamic Flows," 

1978 lecture series, 4, VKI for Fluid Dynamics, March 1978. 

[6] Nodin, F. , Projet de fin d'fitudes, Sup. A6ro [Final Degree 
Project for the Ecole Sup^rieure de 1 ' Afironautique] , 1979. 

[7] Werle, H., "Ecoulements d6coll6s--Etude ph6nom6nologique h 
partir de visualisations hydrodynamiques" [Separated Flow; 
Phenomenological Study Based on Flow Visualizations] , AGARD 
Conference Proceedings No. 168, (1975). 

[8] Legendre, R. , "La condition de Joukowski en Scoulement tridimen- 
sionnel" [The Joukowski Condition in Three-Dimensional Flow] , 

La Recherche Aerospatiale , 5 (1972) , 241-48. 

[9] Durand, E., Electrostatique et magn6tostatique [Electrostatics 
and Magnetostatics] , Masson et Cie, 1953. 

[10] Jacob, C. , Introduction mathdmatique A la M6canique des Fluides 
[Mathematical Introduction to Fluid Mechanics] , Bucharest, 1959. 

[11] Milne-Thomson, L.M., Theoretical Hydrodynamics , 3rd edition, 
Macmillan, 1952. 

[12] Sarpkaya, T. and C.J. Garrison, "Vortex Formation and Resistance 
in Unsteady Flow," Transactions of the ASME- Journal of Applied 
Mechanics (March 1963) . 

[13] Carriere, P., "Notions gSndrales sur les problfemes de confluence" 
[Some General Comments on Confluence Problems] , La Revue Franyaise 
de Mgeanique , 24 (1967) . 

[14] Billet, Thfese de doctorat, a paraitre [Forthcoming doctoral thesis] 

[15] Chometon, F. and P. Fontanet, "Etude thSorique et expdrimentale 
de I'Scoulement bidimensionnel ddcolld autour d'un profil de 
v^hicule automobile" [Theroetical and Experimental Study of Two- 
Dimensional Flow Past the Profile of an Automobile] , 13th AAAF 



Colloquium, Lyons, 1976. 

{16] Siriex, M., "DScollement turbulent en ^couleraent bidimensionnel" 
[Turbulent Separation in Two-Dimensional Fluid Flow] , AGARD 
Conference Proceedings No, 168, (1975). 


32 




